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SUMMARY 


Results are given of an analytical and experimental investigation 
of stresses in .symmetrical .Butt Joints. The methods of analysis, • 
which are Based on the preliminary investigation of NACA TN No. 1051, 
are presented in the first part of the paper. A recurrence formula 
is developed which in conjunction with the appropriate Boundary 
equations can Be used to obtain sets of simultaneous linear equations 
the solutions of which result in the Bolt-load distribution for 
joints of uniform dimensions with Bolts, in line with the load. A 
procedure is also given in which the recurrence formula is applied 
as a homogeneous finite difference equation of the second order. 

In addition, an approximate analysis Based on the shear -lag solution 
of a substitute single stringer structure is presented which may Be 
employed in most practical designs with some gain in time at a small 
sacrifice in accuracy. An example is solved to demonstrate the use 
of the shear -lag solution and a comparison is made with the other 
methods of analysis. 

The second part of the paper describes strain-gage tests for 
joints with five and nine Bolts in line. (The preliminary investi- 
gation analyzed joints with only three Bolts in line.) Because of 
the generally satisfactory agreement obtained in these static tests , 
it appears probable that this analysis may serve as an adequate 
Basis of design of joints .subject to fatigue loads. 


INTRODUCTION 


In a preliminary investigation of Bolted joints the inadequacy 
of the elementary engineering formulas for the stress analysis of 
Bolted joints within the elastic and plastic ranges. But excluding 
failure, was clearly indicated. In this investigation (reference 1) 
a method was developed for calculating the loads carried By the 
individual Bolts in symmetrical Butt joints. The general Bolt -load 


2 


NACA TN Ho . 1^58 


behavior in the elastic range for joints was given as an equation 
which expresses the relationship between the loads on any two 
successive bolts in line with the load. Although this expression 
is applicable to joints of variable bolt spacings, stiffnesses, and 
materials, its application is somewhat tedious, especially in the 
case of long joints. 

The present paper contains the development of a recurrence 
formula for the bolt loads for the simplified case of a symmetrical 
butt joint with bolts spaced evenly in line with the applied load. 

The method presented herein is based on the fundamental relationship 
which was developed in reference 1 between the loads on any two 
successive bolts. The recurrence formula together with the 
appropriate boundary equations furnishes- the stress analyst a 
simple method for the basic analysis of joints . 

In order to reduce the amount of computation .involved in the 
stress analysis of relatively long joints of uniform dimensions, 
the recurrence formula is also treated as a second order finite 
difference equation with constant coefficients. Application of the 
solution of this equation results in a simple and direct determination 
of the bolt loads. This solution, of course, .may also be readily 
applied, to short joints of constant parameters. 

In addition to these methods of analysis an approximate procedure 
is developed based on the shear-lag theory pf reference 2 . A 
"substitute joint", such as that used in reference 3 with a continuous 
connection between plate and straps instead of connections at discrete 
points is analyzed. Equations analogou.3 to those U 3 ed in shear— lag 
problems are derived and a numerical comparison is made with the 
solution of the finite— difference equation. 

The present paper also gives the results of an experimental 
investigation conducted to substantiate further the adequacy of the 
elastic theory as well as to yield additional data on the critical 
bolt load and the behavior of long joints in the plastic range and 
at the ultimate load. The critical load as used in the present paper 
is defined as that bolt load at which either yielding of the plates 
in bearing under the most heavily loaded bolts or yielding of the 
bolts In shear or bending occurs. The test specimens were symmetrical 
five— and nine— bolt joints. A summary of the test data. from the 
present investigation and those of reference 1 is made to help 
formulate principles for design above the limit of elastic action. 
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SYMBOLS 


A cross-sectional area, square inches 

C bolt constant, dependent upon elastic properties, geometric 

shape, dimensions and manner of loading of bolts, and 
upon bearing properties and thickness of plates, incnes per 
kip 

D bolt diaraeter, inches 

E Young’s modulus, ksi 

G shearing modulus of elasticity, kei 

I geometric moment of inertia of bolt, inches^ 1 

K plate constant for tension or compression loading, dependent 

upon geometric shape, dimensions, elastic properties of 
plates, and assumed stress distribution, inches per kip 

L length of joint, inches 

P external applied load, kips 

R bolt load, kips 

b plate width, inches 

k shear— lag constant, dependent on pitch. Young* s modulus, 

bolt constant, and plate areas 

p pitch, inches 

q shear flow, kips per inch 

t thickness, inches 

x distance measured along longitudinal axis of joint 

a, 3 coefficients in finite-difference equation 

& total longitudinal displacement between main plate and butt 

strap, inches 

A plate deformation 


b 
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X = cosh 


tensile strain 
-1 


0*0 


rp = 


2Kp + K s 


t shearing stress^ ksi 

Subscripts: 
av average 

b bolt 

bb bending of bolt 

br bearing 

bs shear of bolt 

cr critical 

1 designation for any bolt 

n number of bolts in joint 

p designation for main plate 

s butt strap 

exp experimental 

theor theoretical 

ult ultimate load 


METHODS OF ANALYSIS 


Development of Recurrence Formula 


Ba sic assumptions of present theory ,— The distribution of loads 
in a bolted joint i3 a statically indeterminate structural problem. 
In order to 3 olve this problem certain basic assumptions and 
definitions must bo made. The assumptions and definitions used 
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herein are the same as those used in reference 1 and are summarized 
as follows: 

(1) The joint is a symmetrical "butt joint where the 
hutt straps are of the same thickness and material. 

(2) The ratio of stress to strain is constant. 

( 3 ) The stress is uniformly distributed over the cross- 
sections of the main plate and "butt straps. 

(4) The effect of friction is negligible. 

(5) The bolts fit the holes initially, and the material 
in the immediate vicinity of the holes is not damaged or 
stressed in making the holes or inserting the bolts. 

(6) The relationship between bolt deflection and bolt 
load is linear in the elastic range. 

General relationship between the loads in successive bolts .- 
On the basis of the aforementioned assumptions, reference 1 shows 
that, for symmetrical butt joints, the general relationship between 
the loads on any two successive bolts in a line with the applied 
load is 



where R is the bolt load; P, the joint load; C.| and C^ + p, 
the bolt constants for the i and i+1 bolts, respectively; Kp 
and Kg are the plate constants for the part of the main plate 


and butt straps between these bolts, respectively; and / R is the 

<r 

sum of the bolt loads Rq to Rj,(See fig. 1 for bolt and space 
designations.) The general expression for the bolt constant C as 
derived in reference 1 is 

r 2t s + ' b p 8t 0 3 + l6t s 2 tp + 8t g tp 2 + t p 3 

3^1) + 192E^bIb 


2t + t 

—2 2 + 

v P \ r V 


®br 




Jjx b: 


( 2 ) 
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and for the plate constant is 


K = 

btE 


( 3 ) 


This analysis was based upon the laws of statics and upon the 
principle of elastic continuity, which requires that after load is 
applied the deflection of holt, i plus the elastic deformation in 
the butt straps between the bolts must equal the deflection of 
bolt i+1 plus the elastic deformation in the main plate between 
the bolts. Equation (l) may be generalised tc apply to joints with 
tapered straps and with bolts of variable spacing and stiffness. 

Derivation of recurrence formula, for joints of constant 
parameters A case that frequently occurs in design is that in 
which all the bolts are of the same material and size and are spaced 
uniformly in line with the applied load. Then 


and equation (1) becomes 



w 


Similarly for bolts i and i-1 



( 5 ) 


Subtracting equation ( 5 ) from equation (h) yields the basic 
recurrence formula for the bolted- joint problem 


‘i-1 





( 6 ) 


For joints with a butt-strap thickness of one -half the main 



2Kp = K s 
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and the recurrence formula "becomes 



tp 

When the joints are made of 24S-T plates with t s = fastened "by 

alloy-steel "bolts, the expression for the "bolt constant (equation (2)) 
reduces to 

C = 3 

Equations ( 7 ) and (8) are directly applicable to the joints tested 
in the present investigation. For other symmetrical "butt-joint 
arrangements, expressions for C similar to equation (8) and "based 
on equation (2) may "be found in appendix A of reference 1. A 
recurrence formula similar to the one in equation ( 7 ) can easily "be 
obtained for joints with varying bolt and plate constants by the use 
of equation (l) directly instead of the simplified bolt-load 
relationships of equations (4) and ( 5 ) . 



Boundary conditi ons . - Before the system of simultaneous equations 
can be solved, the boundary conditions at the ends of the plate must 
be defined. In the joint shown in figure 1(a) the applied load is 
distributed through a fork-like fitting which consists of a main 
plate and two butt straps'. The boundary equation for the left end of 
the joint is 




Rq + El 



(9a) 


and for the right end is 




(9b) 


These equations were derived in a manner similar to the general 
bolt-load relationship in appendix A of reference 1. 


The joint shown in figure l(b) is composed of two idential 
plates carrying equal loads and separated by a filler or idler 
plate connected to the main plate by a row of bolts. For this case 
only one boundary equation is required owing to the double symmetry* 
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of the joint and is given as 


1 + 


(^) 



( 10 ) 


Solution of recurrence formulas and boundary condition equations . - 
By the use of the recurrence formula given in equation (7) and the 
"boundary equations (9) , a system of n simultaneous linear equations 
involving n unknown "bolt loads is obtained. A rapid and accurate 
solution of these equations may he obtained by using the Crout method 
described in reference 4 . When the joint parameters are variable, 
however, the system of simultaneous equations resulting from the 
application of equation (l) can be solved more rapidly by the use 
of the method presented in reference 5. This method takes advantage 
of the analogy between these simultaneous equations and those 
obtained for the current distribution in a direct -current network. 


Since the recurrence formula previously shorn in equation (6) 
is a hemogeneous finite -difference equation of the second order with 
constant coefficients, a solution to this equation may be obtained as 
described in reference 6. Application of the solution results in a 
veiy simple and direct determination of the bolt-load distribution 
in joints of uniform dimensions. The solution of equation (6) is 
shown as 


Solution of Problem by Means of 


Finite -Difference Equation 


- ae U + 3e' U 


( 11 ) 


where the exponent X may be obtained from the relation 



( 12 ) 


where 



( 13 ) 
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The constant coefficients a and 3 of equation (11) are 
determined by the use of the boundary equations (9), as shown in 
appendix A. 

The results are 


a = — 


p = 


where n denotes the number of bolts in the joint. With the 
constants a, p, and X determined., the load carried by each bolt 
my be found directly by successive application of equation (11). 

In appendix B, a numerical example using this method of analysis 
is given. 


Approximate Solution by Shear— Lag Analogy 

C cmparlGon between th e bolted— join t proble m a nd the shea r— lap: 
•problem. — The fundamental action of a bolted joint under load closely 
parallels the action in skin and single stringer combinations used 
as axially loaded panels. (See reference 2.) In both cases the axial 
load is transferred from one component of the structure to another 
through a shear carrying medium. The difference between the two 
actions is that in a bolted joint the loads are transferred in finite 
amounts at definite points rather than through infinitesimal elements 
a3 in the single— stringer struc ture . In che bolted joint, moreover, 
the defoimations of the connecting agent are not solely due to shear, 
but the bolts, being discrete connectors, deform by bending and 
bearing as well as shearing action. Therefore, in order to apply the 
basic equations of the shear— lag analysis, a "substitute structure” 
must be used. 

Expressions for bolt— lead distr ibution base d on modificat ions 
of the shear— lag theory .- An actual joint (fig. 2(a)) may be idealized 
as shown in figure 2(b). The substitute struc ture is obtained by 
distributing the bolts, or shear-carrying medium, over the pitch 
distance p as a "cementing layer" and considering the resisting 
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shear flow q in this material to he 



This expression is analogovis to the shear flow 'rt in the shear-lag 
problem. AI30 , in conformance with the basic assumptions previously 
outline!, the elastic deformations of this fictitious continuous 
cementing medium are assumed proportional to the bolt load; that is. 


where 6 is now defined as the total longitudinal displacement 
between the main plate and strap. Equation ( 15 ) is analogous to the 
shear-strain relationship (7b) given in reference 2. 

By substituting these "equivalent" expressions into the funda- 
mental equations of the shear-lag analysis, the equations used for 
the solution of the bolt— load distribution in bolted joints of 
constant dimensions are obtained. For the practical case shown in 
figure 2 the following equation applies 


shear— lag parameter appearing in analytical solutions for single- 
stringer structures, is defined by 


When the butt-strap thickness equals one-half the main plate 


5 = CR 


( 15 ) 


2 




( 17 ) 
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and equation (l6) reduces to 


cosh It 




1 -^ 


2 sinh 


KL 


(18) 


where t is now defined by 

k 2 - -~ 2 - 
pECA s 


( 19 ) 


It is convenient to take the positive ;c-direction as shown in 
figure 2, starting at a distance of one-half the pitch from the first 
holt. Thus, the length of the joint L may he considered simply to he 
equal to the number of holts in the joint times the pitch (L = np). 


These expressions , strictly speaking, are only "exact" for joints 
fastened hy holts spaced infinitely close together. The accuracy of 
this method when applied to joints with pitches of finite length, 
however, can he shown hy a numerical comparison with the exact 
solution of the finite-difference equation. Appendix 3 illustrates 
the application of the method; and a comparison for a nine-holt joint 
is made of the results obtained hy use of the three solutions presented 
herein. For this case the holt loads determined hy using the shear- 
lag analysis are only about 2 percent less than those computed hy the 
exact solution. 


EXPERIMENTAL INVESTIGATION 
Test Specimens and Procedures 


Specimens .— Tests reported in reference 1 had been confined to 
short joints with two and three holts with a large pitch. In order 
to obtain experimental data on longer joints with a smaller pitch and 
a greater plate— thickness range, tests were performed on six 
symmetrical butt-joint specimens. 


The specimens were constructed of 24S— T aluminum-alloy plates 
fastened by ~ -inch aircraft holts. All specimens were made symmetrical 


about their longitudinal center lines. Only two, however, were short 
enough to he tested as doubly symmetrical joints, that is, symmetrical 
also about the transverse center line. In all joints a pitch 
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of inches was used. Tills pitch was determined by the minimum 

space needed to accommodate the strain gages. The butt— strap 
thickness in all spec imens was one— half of the thickness of the 
main plate. The six specimens were divided into two groups of 
three” joints each. The joints of one group, group C, had five holts 
and those of the other group, group D. had nine holts. (Groups A 
and B were those of reference 1, which are included again in this 

paper. ) The specimens of group C had a width of lg inc nes ^ whereas 

those of group D had a width of 3^ inches. A tabulation ot the 

dimensions of the specimens of groups C and D is shown in table 1 
and photographs of the fractured specimens are shown as figures 3 
and 4. 

In each group of specimens there was one joint for each of the 
three cases found in actual structures. Specimens C--1 and D—l were 
of balanced design, specimens 0—2 and D--2 were designed so that the 
bolts would fail in shear, and specimens 0-3 and D-3 were designed so 
that the plate would fail in tension. All the designs were based on 
the usual assumption that the load is divided equally among the oolbs. 
The same precautions that were taken in the investigation of refer- 
ence 1 to eliminate bearing of the plates on the bolt threads and 
friction of the nuts on the plates were observed. 

Testing procedure.- The test setup of a typical specimen is 
shown" in figure" 5. The joints were tested in tension by means of a 
hydraulic testing machine having a 300-kip capacity and an accuracy 
to about l/2 percent. Strains' were measured on the butt straps with 

electrical resistance- -type gages of f -inch gage length in 1 6 to 20 

increments until failure occurred. Two gage patterns were used as 
shown in figure 6, the second pattern having been considered more 
suitable for the joints with the wider plates. Each specimen was 
preloaded three times to approximately 50 percent of the estimated 
ultimate load. 

Calculatio n of bolt, loads from strain data.- The load on any 
bolt was considered to be the difference between the loads in the 
butt straps at sections midway between the bolt in question and it's 
two adjacent bolts. In order to stud 1 / the influence on the bolt loads 
of the method used to obtain the but b-s trap loads, two independent 
methods were used to convert the strain- -gage data to butt— strap loads. 

By the first method, the butt-strap load was computed simply 
as the gross area of the butt strap times the average stress . The 
average stress was considered to be the product of 'the arithmetical 


NACA TN Wo. IU58 


13 


average of the five strain readings on a gage line and the modulus 
of elasticity (assumed to he 10,600 ksi). By the second method the 
hutt-strap load was computed hy multiplying the area under a curve 
formed hy connecting the five strain readings with straight lines 
hy the thickness of the strap and the modulus of elasticity. This 
area was found hy using the trapezoidal rule. 

Since the trapezoidal method approximates an integration of 
strain across the hutt strap, it naturally is the more accurate 
method. The loads computed on the basis of average strain, however, 
corresponded closely to the one s computed hy the trapezoidal method 
except in some instances where the variation of strain across the 
cross section of the strap was large. In all cases the strain 63 

measured directly in line with holts tended to read lower than the 
outside strain. This tendency was accentuated at higher loads when 
the holt began to hear against the plates, this effect results for 
some cases in a change of strain from tension to compression. Even 
with these large variations in strains the greatest difference in loads 
oalaulated hy the two methods was about 23 percent and this difference 
occurred at a critical load. At lower Joint loads all differences 
were smaller . 

It is also of importance to note that strain measurements taken 
at the center line of the Joint and reduced to load hy the trapezoidal 
rule and compared to the machine load indicated that the internal 
load in the strap was determined within about 5 percent of the 
actual Joint load. Curves representing this relation between' the 
applied Joint load and the measured internal load were linear up to 
Joint failure . This linearity proves that the presence of the lateral 
bending of the hutt strap due to eccentric loading that was evident 
in joints tested in the preliminary investigation was entirely absent 
or negligible in the present tests. The elimination of bending in 
these joints may he attributed, to a large extent, to the fact that 
a greater number of holts were used to resist load and also that the 
increased plate width of these specimens afforded greater flexural 
resistance. Since no correction of the plate loads is necessary, 
these plots are not shown. 

The plotted points in figures 7 to 15 are based upon the 
trapezoidal rule. The analytic curves shown in conjunction with 
these points, however, assume the stress to he uniformly distributed. 
To modify the theory for the irregularity of stress caused hy stress 
concentrations, large hearing deformations and effect of small 
pitches would involve a correction of the plate constants Kp and 
which would result in nonlinear curves of joint load against holt 
load for all holts of a joint. An approximate analysis assuming -the 
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ctress distribution found in specimen D— 1 indicates that, in the 
holt carrying the greatest load, a 14 percent Increase in K increases 
the bolt load only 5 percent. It is apparent, therefore* that the 
calculations for the plate loads are not sensitive to small changes 
in the plate constants. For the other Joints, calculations made by 
use of stress patterns typical of each specimen indicate thi3 same 
tendency to an even greater extent. 


Elastic Behavior 

Curves showing the relationships between the Joint load P 
and the bolt loads B for all test Joints are shown in figures 7 
to 12. Plotted for comparison are the analytical curves, which are 
shown only up to the load above which they are no longer considered 
applicable. The calculated bolt and plate constants and analytical 
bolt loads based on measured dimensions are shown in table 2. 

Figures 13 to Ip show the theoretical and experimental, bolt— load 
distribution for each Joint at the load at which the critical bolt 
load E „ was reached. The bolt load is expressed as a dimensionless 

ratio of the bolt load to the average bolt load P/n. 

The agreement between the experimental points and theoretical 
values can easily be seen in figures 13 to Ip. The general trend 
of the experimental test points follows the theoretical curves; 
however, in some of the specimens there arc discrepancies in the 
individual test points as high as pO percent. In most of the canes 
where these large errors appear, an adjacent bolt has an error of 
approximately the same amount but of the opposite sign. This 
interchange of load is presumably due in large part to irregularities 
of fabrication. If a bolt does not fit tightly it will not "pick up" 
its share of load and the lead that it docs not pick up will, be 
taken by the adjacent bolts. 

On specimens D-l and D-2 bake 11 te gages were used. These gages 
did not adhere well, and. the results obtained with them are open to 
considerable doubt. On specimen D-2 the gagos became detached 
completely during the test, and the last five bolt loads could not 
be determined. 


Inelastic Behavior 

Determination of the "crit i cal bolt load ” Examination of the 
data obtained from the present tests shows that, as was the case in 
the previous investigation, there is seme definite load for each 
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joint "beyond which elastic action no longer continues . The yielding 
of any component , either plate or bolt, is considered to constitute 
the beginning of the inelastic action of the whole joint. Yielding 
in 1 regions of stress concentrations, however, is not included 
as such yielding has no appreciable effect on the over-all elastic 
behavior of the joint. 

The sane method that was described in reference 1 is used 
herein to obtain the limit of elastic action known as the critical 
bolt load. Examination of the test data shown in figures 7 to 12 
shows that the test points for all bolt-load curves tend to break 
away from the linearity of the lower part. The bolt load at the 
intersection of the straight-line portion of the lower part of the 
P-E curve with that of the upper part determines the critical bolt 
load. The curve of the bolt carrying the greatest load when yielding 
occurs is always used to determine the critical bolt load of a joint. 

In table 3 the critical loads and stresses for tests of the 
present investigation are given and in table 4 all the critical loads 
for the tests of this series and also for those of reference 1 are 

listed with their corresponding and g-ratios. A comparison 

of the critical. -bolt-load values shows that although it appears 
generally true that the critical bolt load is largely dependent on the 
parameter D/ 1 0 (as was indicated by Volkerson in reference 3 ) it 

also is influenced to some extent by the ratio of b/D. 

Behavior of bolts for loads above R cr and at failure.- In the 

preliminary work on two- and three -bolt joints (reference 1), it was 
observed that, for loads above the critical bolt load a process 

of bolt -load equalization took place as a result of the yielding of 
plates and bolts so that, at failure of the joint, the bolts carried 
equal loads. This phenomenon was represented by a straight line 
connecting the point representing Per with the point plotted f or 

the average bolt load at joint failure. Examination of the test 
results shows that the same general tendency for the bolt -load curves 
to depart from linearity at R cr that was found in reference 1 is 
seen in plots of the present tests. The bolt load3. however, do not, 


of the bolts to equalize their loads can be explained in the same 
manner as was the failure of the experimental bolt loads below R cr 
to agree with the theoretical values. Examination of the plots in 
figures 7 to 12 shows that in most cases the bolt loads that did not 


in general, approach the average load at failure 



This failure 
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agree -with the theoretical values in the elastic range were also in 
disagreement with the average load close to the ultimate values. 

Despite the fact that the loads in the individual bolts did 
not approach P/n at loads Just below failure, satisfactory 
agreement was found between the observed ultimate and the calculated 
ultimate joint loads baaed on the conventional method of design 
except for the joints of balanced design. In table 4 a comparison 
between the observed and calculated ultimate loads is made for all 
groups. In the calculation for the ultimate loads an allowable 
shear stress of 83 ksi was used for the bolts. This allowable 
shear stress was based on the failing stresses of nine aircraft 
bolts. In order to include the effects of stress concentrations 
and filled holes, the allowable tensile stresses were determined 
from a number of riveted joints with different ratios of b/D. 

The allowable stresses used are based on the ultimate tensile stress 
of standard tensile specimens of 24 S— T .with solid cross sections 
(70 ksi) corrected for these effects. For these joints, with ratios 
of b/D of 7 • 5 , ana 14 , the allowable stresses were taken 
as 66.7, 65. 3 , and 60.7 ksi, respectively. In addition, a value of 
90 ksi was used for the bearing allowable stress, as stipulated in 
reference 7. 


CONCLUSIONS ' 


The analysis of the data from tests of twelve symmetrical butt 
joints (including six of NACA TN No. IO5I) ma.de of 2 kG-T aluminum- 
alloy plates joined by either two, three, five, or nine alloy-steel 
bolts of- the same size with the bolts in line with the axial load, 
leads to the following conclusions: 


1 . The analytical formulas presented were adequate for describing 
the action of these joints in the elastic range because, in general, 
the differences between the test results and the calculated results 
foi" the maximum bolt loads are smaller than scatter of test results 
caused- by uncont: reliable irregularity in the behavior of the 
structures. 


2. The ultimate strengths of the test joints 

(ratio of bolt diameter to plate thickness — - - 

*P 


with thin plates 
1.34 to 3.12) were 


predicted within about 4 percent by the usual assumption that the load 
is uniformly distributed among the bolts. For joints with thick plates 
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ft' 


*) 


- 0.33 "to O .50 ) the prediction, "based on the assumption was 
about 3 percent conservative. 


3- I 11 the ultimate strength calculations of the balanced -design 


Joints of this investigation 



= O.67 to O'. 80' 




however, the 


prediction was about 1 ? percent unconservative, this result 
indicates that the question of determination of the failing loads is 
by no means settled. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., July 2 9 , 1947 
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APB5EDIX A 

DETE3MTI'IA.TION OF COEFFICIENT?.? USED IN 
FINITE— DI]?FEEENC31-EQUATI ON SOLUTION 

When the general recurrence formula (equation (4)) is considered 
as a homogeneous finite-difference equation of the second order,, the 
solution is 

% - oe Xi + 0e -;vi (Al) 


where 


and 



q> = 


2KL + K 

jj b 


For the butt joint, the expression for the left boundary condition 
may be given in the form of equation (9a) as 

2K 0 

-(1 + (?)E j + ^ = - — P (A2) 

When equations (Al) and (A2) (i = 0 and i - 1) are combined, the 
result is 


2KL 


(1 + cp - e^)a + (1 + <p - e“^)p - ~ P (A3) 


The equation expressing the condition that the summation of 
the internal bolt loads must equal the applied load is simply 


n-1 

0 
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or, in expanded form. 


r - 

) - - 

a + - 

(e- n ^ - l) 

(x 

( —X 

\ e - V 


(e - 1) 


(Alt) 


Solving equations (A 3 ) and (A4) simultaneously gives the arbitrary 
coefficients a and ( 3 , Thus, 


a = 



P 


and 
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.APPENDIX B 

NUMERICAL EXAMPLE OF METHODS OF ANALYSIS 


As a numerical example to illustrate the application of the 
three methods of analysis developed, the analysis of test specimen 
D-l will be given. This nine-bolt joint is made up of the following 
components : 

Steel bolts: 

1 

D = -r in. 

4 

= 29,000 ksi 

24S--T plates: 

t = .? in 

JP <3 

t_ = in. 

3 16 

P = l| in. 

* - 3 t ln - 

E = 10,500 ksi 


Analysis by Recurrence Formula 


thickness 


Since the butt— strap thickness equals one— half the main plate 

/ t^\ 

in this joint, the recurrence formula 


V*8 = f 


(equation (6)) applies: 

( ** 

% -l - (? + + B 1+1 - 0 


(Bl) 
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The plate constant is 

K s = -2- 1^25 „ _1_ 

tt s E (3. 5) (0.1875) (10500) 5514 


The holt constant may he determined from equation (8)5 thus 



and 


C 

K_ 


C 


SIC 


0 


1 

433 


433 

5514 


0.07855 


0.1571 


With these coefficients determined, the system of equations 
found according to equation (Bl) and the appropriate boundary 
equation (9b) is 


- 1.1571 K 0 + E-, = -0.07855 P 
E 0 - 2.1571 E 1 + P 2 = 0 
E x - 2.1571 E 2 + E 3 = 0 

e 2 - 2.1571 e 3 + E^ = 0 

E 3 - 2.1571 R h + E^ = 0 
e 4 - 2.1571 e 5 + Eg = 0 

E^ - 2.1571 Eg + By = 0 

Eg - 2.1571 By + Bo = 0 

Ey - 1.1571 Eg = -O.O7855 P 
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The solution of this system of simultaneous equations yields the 
holt loads carried by the individual bolt3. Inasmuch as this set 
of equations is symmetrical about the middle equation, only the 
first five expi’essions need be used. The results are listed in 
table 5* 


Comparing the bolt loads computed by this procedure with the 
bolt load predicted by the conventional analysis in which each bolt 

is assumed to carry the same load (n = — ) shows that the two end 


bolts are overloaded and that the interior bolts carry less than 
they are considered to support. Thus, 


B " B 


= 1.57 


3 

B B 


,11 


R 


2 

T> 

E 

E, 


= — = 0.83 

E 

E 5 

= 0.68 

E 


= 0.63 


Analysis by Solution of Finite-Difference Equation 

Equation (ll) i3 the closed-form solution of the recurrence 
formula applied as a finite-difference equation and is given in the 
form 

Ej^ = ae^ + |3e (B2) 


Since for this case. 


f? 5? 

c " = c 


0.07855 



= 0.1571 
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therefore, from equation (12) 


X = cosh 1 ^1 + 2^ e cosh" 1 ( I.O7855) = O.394 

The coefficients a and ft are determined from expressions (12) 
and (13) , respectively . Since n = 9 for thi3 case, the results are 


a = 


2K 


<p + _*( e nX - l) 


,(e n; ' 

O.157I + O.O786 (e~ - l) 

p^7“3.55)( e ^.39~3 


P = O.OO713 P 


P = 


/ n>. \ 

cp + ----- (e — l) 

_(? r - e“ n 'V (e X - l) 


0.1571 + O.O786 (0^*55 


P = O.1677 P 


Starting with the first holt, successive expressions for each unknown 
holt load are written in the -following manner hy means of equation (32) 

E 0 = Eg = a + p = O.CO713P + O.1677? = 0.1748P 

E 2 = E ? = ae x + pe" X = O.OIO58P + 0.1131P = O.1237P 

E 2 = Eg = ae 2X + pe” PX = O.OI569P + O.O762GP = O.O92OP 

E„ = E = ae 3X + pe~ 3X = O.O2323P + O.O5155P = 0.0748P 

1 5 | 

E^ = ae 4X + pe” 1 ' 1 = 0.03468P + O.O3A60P = 0.0694P 
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By this analysis the individual holt loads are determined directly 
from a 3 ingle expression, obviating the solution of simultaneous 
equations. The results are tabulated in table 5 for comparison 
with the other two methods. 


Analysis by Shear-Lag Analogy 

•Since in this joint A s = A^, equation (lo) may be used 
in the analysis by shear-lag analogy. Thus, 


cosh lc^~ — x) 
R> = kp P 


2 sinh 


kL 


(B3) 


According to equation ( 19 ) j the modified shear— lag parameter k 
is 


r. 



2 2 X 433 

pECA s (1.25)(10500)(3.5)(0.1G75) 


0.1004 


or 

k = 0.318 

The positive x-direction is taken as shown in figure 2; therefore, 
the length of joint L may be considered as 


L = np = (9) (1.25) = 11.25 

By applying equation (B3) successively, the expressions for each 
bolt are given a.s follows: 


E, 


0 


Ec 


cosh O.318 

(o.3i8)(i.25) (o.ViSHii.af) p ’ oa748p 

2 sinh x 


E i = E 7 
E 2 = E 6 


0.0684 cosh 0.318(6.675 - 1.5 X 1.25)P = 0.1230P 
0.0684 cosh 0.318(6.675 - 2.5 x 1.25)P = O.0916P 
0.0684 cosh 0.318(6.675 - 3.5 x 1.25)P = 0.0742P 
0.0684 cosh 0.318(6.675 - 4.5 x 1.25)P = 0.0684P 
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TABLE 1 

co 

CT> 


ELEMENTS OF TEST JOINTS 

2 ij 3 — T plates; S.A.E. 233O (or. equivalent) bolts] 


1 

Number of 

i-in. 

4 

bolts per 
joint 

Nominal 

dimensions 

Measured dimensions 


Specimen 

b 

D 

D 

D 

*8 

t s 

(in.) 

(in.) 

(in.) 

S 

(in.) 

A s 

(sq in. ) 
(a) 

. H 

(sq in.) 
(a) 

Eemarks 

— 

C~1 

{ 


7.5 

0.668 

1.34 

O.183 

O.374 

1.875 

1.875' 

0.344 

. O.7OI 

Balanced design 

C— 2 

5 

7-5 

.334 

.67 

.380 

.749 

I.856 

1.859 

.705 

1.408 

Joint designed to 
fail in bolt shear 

C-3 


7.5 

3.09 

6.25 

.0392 

.0803 

1.874 

1.875 

' .075 

.1506 

Joint- designed to 
fail in tension 

D— 1 


14.0 

.668. 

1.34 

.’89 

.377 

3.501 

3.502 

.662 

1.320 

Balanced design 

D— 2 

9 

14.0 

.334 

.67 

.401 

• 751 

3.502 

3.501 

: 1.406 

' 2.630 

Joint designed to 
fail in bolt shear 

D-3 



14.0 

3.09 

6.25 

.0393 

.0805 

3.5OO 

3.502 

.138 

.282 

Joint designed to 
fail in tension 


®tlross area = b^ 
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TABLE 2 

BOLT ADD PLATE CONSTANTS AND ANALYTICAL BOLT LOADS 
Based on measured dimensions: 


Specimen 

D/tj, 

c 

K s 

S 

Ri/P 

Rg/P 

r 3 /p 

E4/P 


*6^ 

R ? /P 

00 

E 9 /T 1 

C-l 

0.67 

1/433 

I/2885 

I/589O 

0.247 

0.174 

0.152 

O.I76 

0.251 





C -2 

•33 

1/181 

1/5920 

I/II696 

.212 

.194 

.188 

.194 

.212 

— 

— 

— 

— 

C-3 

3-12 

1/151 

1/617 

1/1265 

.270 

.161 

.131 

.163 

.275 





D-l 

.67 

1/433 

1/5581 

1/11090 

.174 

.124 

.092 

.075 

.070 

0.075 

0.092 

0.124 

0.174 

D -2 

•33 

1/169 

I/II806 

1/22088 

.127 

.116 

.108 

.103 

.101 

.102 

.106 

•113 

.124 

D-3 

3.12 

1/152 

1/1156 

1/2368 

.202 

.125 

.082 

.060 

.053 

.060 

.084 

.128 

1 

| .206 
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TABLE 3 


DO 

00 


EXPERIMENTAL LOADS AND STRESSES AT R_ AND FAILURE 

cr 


Specimen 

Joint load 
at critical 
bolt lead 
(kips) 

Critical 
bolt load 
(tips) 

Stress at critical 
bolt load 
(ksi) 

Joint load 
at failure 
(kipe) 

Average stress at 
failure of Joint 

(kBi) 

Type and location 
of failure 



Bearing 

Shear 

Tension 

(a) 

Bearing 

Shear 

Tension 

(*) 


C-l 

21.40 

*5.00 

53.5 

51.0 

35-2 

35-13 

75.3 

71.8 

57.8 

Tension; at bolt 10 
through net section 
of main plate 

C -2 

30.40 

5.20 

24.9 

49.4 

25.O 

42.50 

45.4 

86.8 

35.3 

Shear; all bolts 

C-3 

5-23 

* 1.41 

69.9 

14.3 

40.1 

8.02 

80.3 

16.4 

62.9 

Tension; through net 
section of one butt 
strap at bolt 5, and 
through other butt 
strap at bolt 6 

D-l 

38.OO 

5.20 

55.3 

53.0 

31.0 

63.50 

74.9 

72.0 

51.8 

Tension; at bolt 1 , 
through net section 
of main plate 

D -2 

42.00 

5.25 

28.0 

53-6 

17.2 

75-00 

44.4 

85.0 

30.8 

Shear; all bolts 

D -3 

8.00 

1.33 

66.4 

13.6 

30.6 

15.05 

83.0 

17.1 

59-0 

Tension; at bolt 9 , 
through net section 
of butt straps 


Computed using net area . 

^ Average of maximun "bolt loads In upper and lower Joint. 
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TABIE 4 


SUMMARY OF CALCULATED AND EZPERBENTAL RESULTS 
FOR 2 , 3, 5, AND 9 BOLT JOINTS 


Specimen 

(a) 

D/tp 

0>) 

Vd 

Observed 

P 

ult 

(kips) 

Calculated 

P 

ult 

(kips) 

(c) 

P ult 

Joint load 
at R 

cr 

(kips) 

R 

cr 

(kips) 

Caic. P ult 

A-l 

o.8o 

5 

15.96 

16.30 

0.98 

7.00 

3.88 

A -2 

•50 

5 

16. o4 

16.30 

•99 

8.00 

4.80 

A -3 

I.52 

5 

10.62 

IO.78 

•99 

— 

— 

B-l 

.67 

5 

23.40 

24.50 

.96 

11.40 

4.l6 

B-2 

• 50 

5 

24.20 

24.50 

•99 

13.56 

4.8o 

B-3 

1.34 

5 

12.02 

12.50 

.96 

8.25 

3-24 

C-l 

.67 

7-5 

35-13 

39.60 

.89 

21.40 

5.00 

C-2 

•33 

7-5 

42.50 

40.70 

l.o4 

15.OO 

5.20 

C-3 

3.12 

7-5 

8.02 

8.60 

•93 

5.10 

1.41 

D-l 

.67 

l4.o 

63.50 

73*40 

.87 

38.00 

5.20 

D-2 

•33 

14.0 

75.OO 

73.40 

1.02 

42.00 

5.25 

D-3 

3-12 

14.0 

15.05 

15.95 

•95 

8.00 

1-33 


Specimens A end. B from reference 1. 

^Based on measured dimensions. 

c Calculated ultimate loads are Based on conventional design method. 
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tabu: 5 

COMPARISON OF BOLT LOADS IN SPECIMEN D-l 
AS FRACTION OF TOTAL LOAD 


Bolt; 

Methods of analysis 


Recurrence 

formula 

Finite -difference 
equation 

Shear -lag 
analogy 

0 

0.1748 

0.1748 

0 .1748 

1 

.1237 

.1237 

.1230 

2 

.0920 

.0920 

.0916 

3 

.0748 

.0748 

.0742 

4 

.0694 

.0694 

.0684 

5 

.0748 

.0748 

.0742 

6 

.0920 

.0920 

.0916 

7 

.1237 

.1237 

.1230 

8 

.1748 

.1748 

.1748 


1.0000 

1.0000 

0 .9956 
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la) Symmetrical bull joint. 


P — 


P/2 

P/2 



lb) Filler-plate joint. 


P/2 

P/2 
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Figure I. — Symmetrical butt joints with bolts in line with the 
axial load . 
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(a) Actual joint. 


P 



— ^ P/2 
P/2 


tb) Substitute joint. 






- P/2 

Cementing 

lager 


Cc) Half - structure . 
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Figure 2.-Axiallg loaded butt joint. 
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Figure 4.- Side view of fractured specimens. 
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Figure 5.- Test setup of a typical specimen 



Cb|- 
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Group C 



Group D 
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Figure 6 ~ Strain gage arrangements across straps 
of test specimens in groups C and D. 



o 



Figure 7— Observed relationships between applied joint load and bolt lood for specimen C-l 
and comparison with calculated values. 
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Bolt load , R , kips 

Figure 6. — Observed relationships between applied joint load and bolt load for specimen C-2 
and comparison with calculated values. 
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Applied joint load , P, Kips 


co 



Figure 9.- Observed relationships between applied joint load and bolt load for specimen 
C-3 and comparison with calculated values. 
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Figure 10- Observed relationships between applied joint load and bolt load for 
specimen D-l and comparison with calculated values. 
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Figure I !.- Observed relationships between opplied joint load and bolt load for specimen D-2 

and comparison with calculated values. 


NACA TN No. 1458 


Applied joint load , P, kips 



Figure IE- Observed relationships between applied joint load and bolt load 
for specimen D-3 and comparison with calculated values. 
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Specimen D-l 



Specimen C-l NATI0NAL A0VIS0RY 
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Figure 1 3.- Comparison between the experimental and theoretical 
bolt-load distribution at R cr for joints of balanced -design. 
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Specimen D-2. 


> 

o 

Ol. 

\ 

< 

CO 

CT 

T3 

e 

a 

> 

a 

f*r 


12 


1 1 — 1 — 

p— 

oooool p 

Boll 

Ob 

1 3 4 


1 

3 = I5.( 

D hips 



\ J 






6 


_ci 

cr 


§ 0.4 


o 


•av 


0 

BoltO 


I L 3 
Specimen C-£. 


^theor>/ Rqv 

^exp/ ^qv 
R = R ( 


Figure l4.-CompQri.5on between the experimentoi and 
theoretical bolt -load distribution ot R cr for joints of over-design. 
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Specimen D-3 



Specimen C-3 


Figure 15- Comparison belween the experimental and theoretical 
bolt-load distribution al R cr for joints of under-design. 


